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i We report on a measurement of D° — D° mixing and a search for 

(jj CP violation in the D° — > K + K~ , n + n~ and _ft' ± 7r =F channels. We use 

D°'s coming from D* + decays, so that the flavour of the D° at produc- 
tion is tagged by the charge of the pion that is also emitted. We also 
use an independent set of _D°'s coming directly from the hadronization 
\Q °f the charm quark, but in this case the flavour of the charmed meson 

is not known. We analyze events collected by the BABAR experiment at 
\q the PEP-II asymmetric-energy e + e~ collider, corresponding to an inte- 

qq grated luminosity of 468 fb" 1 . We measure the mixing parameter value 

O to be yep = [0.72 ± 0.18(stat) ± 0.12(syst)]%, and exclude the no-mixing 

hypothesis at 3.3<r significance. We find no evidence of CP violation, ob- 
serving AY = [0.09 ± 0.26(stat) ± 0.06(syst)]% which is consistent with 
:^ zero. 
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1 Introduction 



Mixing in the charm sector is a well-established phenomenon [H El El HI El E] although 
there is no single measurement that exceeds 5a significance. Recently the LHCb [7j 
and CDF [8] Collaborations have reported evidence of CP violation (CPU) in the 
difference of the integrated CP asymmetries in the D° — > K + K~ and D° — > 7t + tt~ 
channels. This result was unexpected at the current experimental precision, and 
it may be a manifestation of New Physics (NP), although a Standard Model (SM) 
explanation cannot be ruled out. These measurements have renewed the interest of 
the community in charm physics as a sector in which to search for NP manifestations. 

Under the hypothesis of CPT conservation the two mass eigenstates (Di and D 2 ) 
can be written in terms of the flavor eigenstates (D° and D°) as: 

\D 1>2 ) =p\D°) ±q\D°) with |p| 2 + |g| 2 = 1. (1) 

If CP\D°) = + \D }, then in the case of no CPU, D x is the CP-even state and D 2 the 
CP-odd state. The parameters that describe D° — D° oscillations are proportional to 
the difference of masses (m^) and widths (Tj) of the mass eigenstates: 

m 1 -m 2 T 1 - T 2 
x = and y = — , (2) 

where T = (Ti + F 2 ) /2 is the average width. Mixing will occur if the mass eigenstates 
differ from the flavour eigenstates, that is, if either x or y is non-zero. SM predic- 
tions for the mixing parameter values are at the order of a percent or less and, at 
present, experimental measurements are in agreement with these predictions. Unfor- 
tunately the theoretical predictions are affected by large computational uncertainties 
on the dominant long-range-diagram contributions, preventing these measurements 
from being strong tests of the SM. 

In the following we present a measurement of the mixing parameter yep [5] and 
the CP-violating parameter AY, defined as: 

r+ + r+ r + - r+ 

y C p=—^ 1 and AY = — — — , (3) 

where T + (T + ) is the average width of the D° (D°) when reconstructed in CP-even 
eigenstates. 

The measured values of yep and AY constrain the parameters that govern mixing 
and CPV in the charm sector. Neglecting the effect of direct CP violation, estimated 
to be at least one order of magnitude below our current sensitivity, we relate yep and 
AF to the mixing and CP-violating parameters as follows: 

Vcp = U cos -^-x sin and AY = — x sin H — ^-y cos 0. (4) 
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The asymmetry Am = (q]p)^+(p/q) 2 measures CP violation in mixing, while <p is sen- 
sitive to CPV in the interference between decays with and without mixing, being 
the weak phase of the quantity A/ = with Af (Aj) the amplitude for the decay 

D°(D°) /. In principle, <j) can depend also on the final state, but with our current 
level of precision we are not sensitive to this 1 1 j . In the absence of CP violation 
Vcp = V and AF = 0. 

2 Data Sample and Backgrounds 

We reconstruct the D° in the h + h~ (h = K, tc) and K ± tt t final states and measure 
three lifetimes: 

• t + for the D° — > h + h~ decays, 

• f + for the D° — > h + h~ decays, 

• for the D° (and D°) — > K ± n T decays (the Cabibbo favored K~ir + and the 
doubly Cabibbo suppressed K + n~ decays are collected in the same sample). 

Due to the small mixing rate (< 1%) we can neglect the effect of mixing and assume 
that all signal proper time distributions are exponential. The untagged D° — > K + K~ 
sample [TTJ is assumed to contain 50% of D° and 50% of D° decays. The three values 
of inverse lifetime are used to compute yep and AF: tk-k is used to access the average 
width T while, r + (f + ) is used to obtain Y + (r + ). 

We use tagged decays of the D° coming from D* + decays, through D* + — > D°ir^, 
as well as untagged decays coming directly from the hadronization of the charm quark. 
The tagged and untagged samples are independent, i.e. sua event containing a tagged 
candidate and at least one untagged candidate is excluded from the untagged sample. 
In the tagged sample the flavour of the D° is determined by the charge of the pion 
that is also emitted. Due to the significantly higher level of background in the 7r + 7r~ 
final state, we do not use the related untagged sample. 

We analyze 468fb _1 of data recorded by the BABAR detector [12] at, and slightly 
below, the T(4S) resonance at the e + e~ asymmetric-energy PEP-II 5-Factory . To 
avoid potential bias, we finalize our data selection criteria, as well as the procedures for 
fitting, extracting statistical limits, and determining systematic uncertainties, prior 
to examining the results. 

An oppositely charged pair of K + or ir + candidates satisfying particle identifica- 
tion criteria is fit to a common vertex to form a D° candidate. We require each D° 
to have momentum in the center-of-mass (CM) frame pqm > 2.5GeV/c in order to 
remove almost completely D 0, s coming from 5-meson decays. For the tagged modes, 
we form the D* + candidate by fitting a D° candidate and a charged pion track 7r+ 
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to a common vertex, which is required to lie within the e + e~ interaction region. The 
7r+ momentum is required to be greater than 0.1 GeV/c in the laboratory frame and 
less than 0.45 GeV/c in the CM frame. We veto any 7r+ candidate that may have 
originated from a reconstructed photon conversion or 7r° Dalitz decay and reject a 
positron that fakes a 7r+ candidate by using energy loss information. We also select 
tagged candidates in a Am window, 0.1447 < Am < 0.1463 GeV/c 2 , where Am is the 
difference between the reconstructed D* + and D° masses. This requirement strongly 
suppresses backgrounds. 

The proper time t and proper time error a t of each D° candidate are determined 
from a combined fit to the D° production and decay vertices. The x 2 probability of the 
vertex fit must satisfy P(x 2 ) > 0.1%. We retain only candidates with — 2 < t < 4ps 
and <y t < 0.5 ps. For tagged decays, this fit does not incorporate any 7r+ information 
in order to ensure that the lifetime resolution models for tagged and untagged signal 
decays are very similar. The most probable value of a t for signal events is ~ 40% of 
the nominal D° lifetime [13] . 

For cases where multiple D* + candidates in an event share one or more tracks and 
the D° decays to the same final state (K~ tc + and K + n~ are considered to be the 
same final state in this context), we retain only the candidate with the highest P(x 2 )- 
If an event contains a tagged decay, all untagged candidates from that event are 
excluded from the final sample. In an event with no D* + candidate and multiple D° 
candidates decaying to the same final state, we retain only the D° candidate with the 
highest P(x 2 )- The fraction of events with multiple D° candidates with overlapping 
daughter tracks is <C 1% for all final states. 

In Fig. [T] we show the reconstructed invariant mass distributions for the selected 
D° candidates in both tagged and untagged modes. We fit the mass distributions in 
order to extract the total number of background candidates. In Fig. [l]we also report 
the fit results and, below each plot, show the normalized Poisson pulls [13]. For the 
tagged CP-even modes, the D° and D° samples are fit simultaneously, sharing all 
parameters except for the expected signal and background candidate yields. 

We perform a mode-dependent, data-driven optimization of the invariant mass 
window position and width, in order to reduce significantly the effects of the linear 
correlation between the reconstructed mass and the reconstructed proper time. The 
signal regions obtained for each mode are shown in Fig. [T] with dashed lines. These 
are 34MeV/c 2 wide for all modes except untagged D° — > K + K~ . In this mode the 
signal region width is reduced to 24 MeV/ c 2 due to the higher level of background, as 
observed in the corresponding plot. We define a lower- and a higher- mass sideband 
each of width 20 MeV/ c 2 ; these are used to study and characterize the combinatorial 
background. The mass sidebands for the untagged modes are ±44.5 MeV/ c 2 from the 
signal region center. In case of the tagged modes, the distance of the sidebands from 
the signal region center is ±35.5 MeV/ c 2 , and the Am window is shifted to higher 
values, 0.151 < Am < 0.159 GeV/c 2 . In the sidebands the tagged and untagged 



3 



■ T 2000 

S 1500 
o 

S 1000 

I 

[S 500 



JIJC 




+2 
-2 

1.82 1.84 1.86 1.88 1.9 1.92 
M„ (GeV/c 2 ) 
xlO 3 



' '.j 2000 

§ 1500 

•J 1000 
'c 

I 500 


g +2 







: D*+ A 




p. rot / \ 









1.82 1.84 1.86 1.88 1.9 1.92 
M„ (GeV/c 2 ) 

xlO 3 




1.82 1.84 1.86 1.88 1.9 
M KK (GeV/c 2 ) 




1.82 1.84 1.86 1.88 1.9 
M KK (GeV/c 2 ) 



1.82 1.84 1.86 1.88 1.9 
M KK (GeV/c 2 ) 






S +2 



; unt 
i Ktz 













1.8 1.82 1.84 1.86 1.88 1.9 1.92 
M K , (GeV/c 2 ) 



1.82 1.84 1.86 1.88 1.9 1.92 
M K , (GeV/c 2 ) 



Figure 1: The reconstructed two-body mass distribution for the seven modes. The 
vertical lines show the signal region. The shaded regions are the background contri- 
butions. The normalized Poisson residuals for each fit are shown under each plot. 



samples are not independent. 

After the selection we divide the backgrounds into two categories. Candidates 
for which the common ancestor of the D° products is a long-living charmed meson 
are collected in the misreconstructed-charm background category. In Table [T] we 
report the composition of this background in the signal region, obtained by studying 
simulated events. The other background candidates, consisting mainly of random 
tracks, fall into the combinatorial background category. In Table [2] we report the 
number of signal and background candidates after selection, for the signal region. 



3 Lifetime Fit 

In order to extract the three lifetime values we perform an extended unbinned- 
maximum-likelihood fit to the 2-dimensional distribution of proper time and proper 
time error. All modes are fit simultaneously: the signal resolution function parame- 
ters are shared among the modes, while the background Probability Density Function 
(PDF) parameters are not. 

The single-mode PDF for the signal events consists of an exponential convolved 
with a resolution function. The latter is the sum of three Gaussian functions with a 
common mean (offset) and widths proportional to the per-event proper time error, 
scaled with three different factors (one for each Gaussian). In order to take into 
account differences in the reconstruction due to the different final states, we multiply 
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57.1 


8.8 


35.8 


a>° ->■ 7t 7t + a:- 


1.1 


70.3 


1.7 


63.3 


6.9 




< 1 


2.9 


< 1 


11.8 


< 2 


D° ->■ K + K~ 


< 1 


< 1 


1.3 


< 1 


3.5 


D° -> 7T + 7T- 


1.8 


< 1 


2.2 


< 1 


3.1 


D° ->■ TT + TT-TT 


< 1 


< 1 


7.0 


< 1 


17.3 


A decays 


< 1 


< 1 


< 1 


4.9 


2.6 



Table 1: Expected composition of the misreconstructed-charm backgrounds. Only 
misreconstructed-charm background modes that have > 1% contribution in at least 
one signal mode are listed. For the tagged modes, the fractions are the sum of the 
separate D° and D° tags. 





Signal 


Combinatorial Bkgd. 


Charm Bkgd. 


Tagged n + n 


65 430 ±260 


3 760 


97 


Tagged K + K~ 


136 870 ±370 


653 


309 


Tagged A^tt^ 


1 487 000 ± 1200 


2 849 


642 


Untagged K + K~ 


496 200 ± 1200 


165 000 ± 1 000 


5 477 


Untagged K ±r n T 


5 825 300 ±2 600 


1 044 552 


4 645 



Table 2: Signal and background yields in the signal region; yields with uncertainties 
are those obtained directly from the lifetime fit to data. For the tagged modes, the 
yields are the sum of the separate D° and D° tags. 

each Gaussian scale factor by another scale factor that depends on the final state 
(the K ± Tt T factor is fixed to 1). In the same way we introduce a third scale factor 
that depends whether the mode is tagged or untagged, fixing to unity the one for 
the untagged modes. For the tagged CP-even modes we also take into account the 
wrongly-tagged signal candidates, fixing the fraction of these events to the value 0.2%, 
obtained from the simulated events. Since the proper time PDF depends on the proper 
time error, we multiply each signal PDF by the 1-dimensional binned distribution 
of a t to avoid biases. The normalization of the proper time PDF is computed for 
each <7 t . The a t histogram for the signal events is obtained from the distribution of 
the events in the signal region after subtraction of the misreconstructed-charm and 
combinatorial-background contributions. 

The 2-dimensional PDF for the misreconstructed-charm background is a signal- 
like PDF, fitted to the simulated events and then fixed in the final fit. Since this is a 
physical background, its lifetime, composition and number of events change with the 
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mass window. Therefore we have decided not to use the sidebands to characterize it. 

The combinatorial PDF is determined as a weighted average of the PDFs in the 
two mass sidebands, which consist of 2-dimensional histograms. For the untagged 
K + K~ mode the sideband PDFs for this category are signal-like. Contributions of 
signal and misreconstructed-charm in the sidebands are parameterized using the sim- 
ulated events, and then fixed. The weighting parameter is determined from simulated 
events and then is varied as part of the systematic studies. 

The expected total-background candidate yields are evaluated from the mass fit 
and then corrected using the simulated-event information. The misreconstructed- 
charm contribution is estimated from the simulated events, and the combinatorial one 
is obtained by subtraction. In the final fit the background yields of the two categories 
are fixed for all modes except for the combinatorial untagged K + K~ mode, where it 
is allowed to float. This became necessary since the prediction of the mass fit was 
not accurate enough for this mode, where the combinatorial background represents 
almost 25% of the events in the signal region. 

4 Analysis Validation and Systematics 

The validation of the procedure has been performed on four independent samples of 
simulated events, each equivalent to data integrated luminosity, and also on a large 
ensemble of pure Toy MC samples. We have also performed a qualitative validation 
on data by running the fit in different configurations. For example, we have fitted 
the tagged and untagged samples separately, finding the K + K~ and K ± n T tagged 
and untagged extracted lifetimes compatible within the statistical uncertainties. We 
let the 7r + 7r~ and K + K~ samples have different lifetimes, allowing for physical effects 
depending on the final states (direct CPV and the dependence of on the final state), 
and found r + and f + to be compatible for the two modes. 

In addition to these tests, we have also identified sources of systematic error and 
have evaluated their contributions, as reported in Table |3j We have evaluated the 
systematic effects due to the choice of the signal region by varying its position and 
width. We have varied the fraction of mistagged events in the D° — > h + h~ tagged 
modes, and the fraction of _D°'s in the untagged K + K~ mode. The proper time error 
PDF is obtained by subtraction of the background distributions. However, in the 
untagged K + K~ mode, the combinatorial yields are extracted from the lifetime fit and 
not known a priori. In the nominal fit we first estimate the number of combinatorial 
events as for the other modes, and use this to perform a first simultaneous fit. We 
then repeat the fit using the yields just extracted, and this fit yields the nominal 
results. In order to evaluate the systematic error associated with this procedure we 
repeat the fit a third time, and take as a systematic error estimate the difference from 
the nominal value. We have varied the misreconstructed-charm lifetimes and yields, 
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charm bkgd. lifetimes 


0.042 


0.001 


comb, yields 


0.043 


0.002 


comb, sideband weights 


0.004 


0.001 


comb. PDF shape 


0.066 


0.000 


cr t selection 


0.052 


0.053 


candidate selection 


0.028 


0.011 


Total 


0.124 


0.058 



Table 3: The yep and AY systematic uncertainty estimates. The total is the sum-in- 
quadrature of the entries in each column. 



estimated using the simulated events, by > 2a. 

The combinatorial PDF is extracted from the sidebands after fixing the signal 
and misreconstructed-charm contributions. We have applied the variations described 
above for the misreconstructed-charm events in the signal region also in the sidebands, 
for both the signal and the charm-background contributions, and re-extracted the 
combinatorial PDF. We have also varied the number of combinatorial-background 
events in the signal region for the modes in which it was fixed, and the weighting 
parameter for each mode. 

We have varied the selection criteria, in particular that on a t by ±0.1 ps. We have 
also estimated the systematic impact of the best candidate selection by removing 
or keeping all the overlapping candidates. We have estimated the effects of SVT 
misalignment and have found these to be negligible. 

5 Results and Conclusions 

The seven projections of the lifetime fit are reported in Fig. [2] The following lifetime 
values are extracted: 

r+ = (405.69 ±1.25) fs, f + = (406.40 ± 1.25) fs, r Kn = (408.97 ±0.24) fs. (5) 

The lifetimes are reported with the statistical error only. The K ± 7r Zf lifetime is 
compatible within one standard deviation with the PDG D° lifetime [H] and the 
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Figure 2: Proper time, t, distribution in each final state with the fit result overlaid. 
The combinatorial distribution (indicated as 'Comb.' in light gray) is stacked on 
top of the misreconstructed-charm distribution (indicated as 'Charm' in dark gray). 
The normalized Poisson pulls for each fit are shown under each plot; "unt" indicates 
the untagged datasets. The bottom-right plot shows the individual lifetime values 
(statistical uncertainty only); the gray band indicates the PDG D° lifetime ±la [13] . 

CP-even lifetimes are significantly lower, as shown in the bottom-right plot of Fig. [2] 
Combining the values of inverse lifetime following Eq. ([3]), we obtain: 

ycp = (0.72 ±0.18 ±0.12)%, and AY = (0.09 ± 0.26 ± 0.09) fs. (6) 

The first error is statistical, obtained from the covariance matrix resulting from the 
fit, and the second error is systematic. This measurement represents the most pre- 
cise measurement of yep, and excludes the no-mixing hypothesis at 3.3a significance. 
The value of yep presented here is compatible with all previous measurements. In 
particular it is compatible with the previous BABAR measurement [3] with a proba- 
bility of > 2%, taking into account that ~ 40% of the events in the current sample 
are also present in the samples used in the previous measurements [21 13] , and that 
the systematic errors are fully correlated. This result favors a lower value for yep, 
and approaches the value of the mixing parameter y when measured directly [15], as 
expected if CP is conserved. We find no evidence of CPV. 

In conclusion, we report evidence of D° — D° mixing with 3.3<r significance, ob- 
taining the mixing parameter value yep = [0.72 ± 0.18(stat) ± 0.12(syst)]%. We 
find no evidence of CP violation, and measure the CP- violating parameter value 
AY = [0.09 ± 0.26(stat) ± 0.06(syst)]%. 
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